Semiclassical Husimi functions for spin systems 



Marcel Novaefl and Marcus A.M. de Aguiar 
Instituto de Fisica "Gleb Wataghin" , Universidade Estadual de Campinas, 13083-970 Campinas-SP, Brazil 

We derive a semiclassical approximation to the Husimi functions of stationary states of spin 
systems. We rederive the Bohr-Sommerfeld quantization for spin by locating the poles of the corre- 
sponding local Green function. The residues correspond to the Husimi functions, which are seen to 
agree very well with exact calculations. 



o 
o 

> 

o 
m 



> 



o ■ 

^ : 

9 L,: 
-i— > ■ 

G ■ 

3 : 



Rigorous derivation of semiclassical approximations in 
phase space via path integrals Q for systems with one 
degree of freedom have recently received considerable at- 
tention, both for continuous variables and spin systems. 
Baranger et al for example, have discussed the canon- 
ical coherent state path integral and its semiclassical ap- 
proximation in some detail, including an initial-value rep- 
resentation and the Green function. The study of semi- 
classical propagation of wave packets, using complex |3| 
or nearly real [4| trajectories, for regular and chaotic 
systems, has developed considerably over the last few 
years. The spin path integral, and its semiclassical ap- 
proximation, has found an important application in the 
study of spin tunnelling and topological effects [g] . Stone 
et al have derived the spin coherent state semiclassical 
propagator in detail pj, paying particular attention to 
the so-called Solari-Kochetov Q correction. This cor- 
rection is related to the difference between the average 
value of the Hamiltonian in coherent states and its Weyl 
mbol 9], and has a counterpart in the canonical case 



To obtain semiclassical approximations for the energy 
levels E n and stationary states (x\n) = ip n {x) of one- 
dimensional bound systems, on the other hand, one nor- 
mally resorts to the usual Bohr-Sommerfeld (BS) and 
WKB theories [ljj. A coherent state version of these 
theories, which works in phase space, is also available 
and produces a BS formula and a semiclassical ap- 
proximation to the Husimi functions H n (z) = |(z|n)| 2 . 
Recently, Garg and Stone have derived a semiclassi- 
cal (BS-like) quantization condition for spin systems, in- 
cluding the first quantum corrections (see also H3)- By 
taking the trace of the semiclassical Green function, they 
obtained the energy levels as the location of its poles. 
In the present work we have obtained the semiclassical 
Husimi functions for spin systems. 

The non-normalized spin coherent states are defined 
by \z) — exp{zJ + }|j, —j), and the semiclassical approxi- 
mation to the propagator K — (zf\e~ lHt ^ h \zi) is 



K sc i(z f ,Zi,t) 



ie iB /i d 2 S \ V2 



h 2j dz t dz f J eXP U^' (1) 



where the phase is the classical action plus an extra term 



'Electronic address: mnovaes@ifi.unicamp. br 



known as the Solari-Kochetov (SK) correction: 



$ = S + I; 



SK 



S 



A(t')dt'. 



The classical spin action is given by 



5' 



ihj 



zz — zz 
l + zz 



H(z,z) 



dt' + B, 



(2) 



(3) 



where the integral is done along the classical trajectory 
determined by the Hamilton equations of motion 



ihz 



1 (III 

g(z,z) dz ' 



ihz = 



1 dH 

g(z,z) dz ' 



(4) 



and the classical Hamiltonian H(z, z) is the average value 
of the quantum Hamiltonian, H(z,~z) — (z\H\z)/(z\z). 
This action obeys the Hamilton- Jacobi relations 



% OS 
h dz~f 



2Mt) 
1 + Zfz(t)' 



i dS _ 2jz(0) ^___ H 
hdzi ~ l + z(0)V ~dt ~ ~ 

(5) 

The function g(z, z), which is unity in the canonical case, 
is given by 



I:) 2 



dzdz 



\n{2 



2j 



(1 + zzf 



(6) 



Note that B = -ihj ln[(l + z f z(t))(l + z(0)zi)] is a 
boundary term, that takes into account the fact that in 
general z is not the complex conjugate of z (the discrete 
time formulation of the path integral indicates that the 
variables z and ~z must be considered as independent, 
so we denote the actual complex conjugate of z by z*). 
That means that if one defines the usual canonical (q,p) 
variables according to 



q + ip 



q-ip 



VI + zz 



sfl + zz 



(7) 



then q and p will in general be complex numbers. 

The semiclassical limit for spin systems consists in let- 
ting h — * and j — > oo, but keeping hj — 1. If the 
Hamiltonian if 0(hj), then S is 0(hj), but the SK cor- 
rection 



A 







1 dH d 



1 



dH 



dz Ag(z, z) dz dz Ag(z, z) dz 



(8) 
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is 0(h), and therefore can be considered small. Note that 
since I^z* the Hamiltonian H(z,~z), the action and the 
SK correction can all be complex. 
The semiclassical Green function, 



G 



scl 



1 f°° 

{z,E) = - K scl (z,z,ty Et / h dt, 
in Jo 



(9) 



can be calculated by making a stationary exponent ap- 
proximation to the integral. Note that we are inter- 
ested only in its diagonal elements. This implies z. t = z, 
~z~f — z*, but in general z(t) ^ z and z(0) ^ z* , so that we 
do not have a real periodic orbit (by real orbit we mean 
one in which q and p are real) . The stationary time to is 
determined by the condition 



d($ + Et) 




dS 




OIsk 


dt 


to 


dt 


H 

to 


dt 



E = 0. 



(10) 



As usual in semiclassical calculations, we do not consider 
derivatives of A, because including such terms would be 
inconsistent with the gaussian approximation involved in 
the derivation of Therefore Ea. I|l()fl can be also writ- 
ten as 



where 



E -£(z,t )+A{z,t Q ) =0, 



£(z,t ) = -- 



(11) 



(12) 



Even with this simplification it is hard to find the poles 
of G(z, E). Garg and Stone |n| have done this by calcu- 
lating its trace under another stationary phase approxi- 
mation, which leads to z(t) — z, "z(O) = z*, and thus to 
real periodic orbits. We take a different route, that will 
allow us to obtain not only the energy levels but also the 
Husimi distributions. Even though the classical orbits 
involved in the calculation of and © are complex, 
we argue that the largest contributions to the function 
G(z, E) (and not only to its trace) must come from the 
vicinity of the real periodic orbit through z, and its rep- 
etitions. The accuracy of the final results support this 
idea. 

Let us denote the period of the orbit through z by T(z) 
and expand the stationary time to as 



t Q ~ nT + T , 



(18) 



where n counts the repetitions of the real periodic orbit 
and To, assumed small, has to be determined. Expanding 
the stationary exponent condition (|10|l we find 



Tn = - 



E-S + A 



where 



:=a(z,nT) = - 

dt 



nT 



d 2 s 

dt 2 



(19) 



(20) 



nT 



is the energy of the classical trajectory, not to be confused 
with E, the argument of the Green function. In order 
to proceed with the integration, we need to expand the 
exponent to second order in time. We define 



d 2 s 

dt 2 



a{z,t ), 



(13) 



to 



and neglect the second derivative of Isk, in order to ob- 
tain 



G(z,E) 



1 

ih 



we tB/hj Q2 S X 

dzidzj J 



1/2 



aj 



exp{-(^}, (14) 



where 



V = S(t )+T SK {t )+Et . 



(15) 



We can find a more convenient way of expressing a in 
order to transform the prefactor. The form (see (J) 



leads to 



G(z,E) = - 
in 



_ d 2 S 

dzfdzi 



1 /^ e ™\ 1/2 



(16) 



zzj 



exp{-(^}. (17) 
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FIG. 1: Exact (top) and semiclassical (bottom) Husimi dis- 
tributions for the simple Hamiltonian H — uihj z . 
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Now we must expand <p to second order in To. Note 
that S{nT) = nS - n£T + B(nT), where 



S = ifoj / ■ = —dt. 

1 + zz 



(21) 



Note also that IgxinT) — 7iTsk(T). It can also be 
shown |2| that 



1 



a( n ) ihg(z, z)\z\ 2 



d 2 S 



(22) 



After n repetitions of a periodic orbit the prefactor ac- 
quires a phase of (—1)". Therefore, the result of this 
expansion is 



tp ~ n [S - £T + ET + T SK - irfi] 



B(nT) - £T 



ET + AT . 



(23) 



If we add and subtract nAT, define x = E — £ + A and 
use T = dS/d£ together with equation l|22|l we obtain 



tp ~ n 



S + Isk — ~kTi 



dS dS, I cPS 2 
d£ X ~ d£ + 2d£* X 



g(z, z) 2ifi\z\ 



B{nT). 



(24) 



We recognize inside the brackets the expression for the 
expansion of S(E + A) around £ . Since A ~ 0(h) we 
further expand 



S(E + A)~S(E) + ^A 
d£ 



(25) 



and we end up with 
(p ~ n[S(E) +1 SK (E) - irk] 
Summing over n we obtain 



/ — N o-Til -12 + B ( nT )- 

g[z, z) 2in\z\ z 

(26) 



G(z,E) 



^{l + \z\ 2 ) 2 3+ l e i[S+IsK-nh]/H 

ihj\z\ l - e i\s+i S K-Kn]/n 
1 (E - £ + A) 2 



x exp 



g{z,z) 2h 2 \z\ 2 



(27) 



where we have used B(nT) = —2ihj\n(l + \z\ 2 ). The 
poles of this function are determined by the condition 



{S + l SK ){E n ) = (2n+l)7r7i, 



(28) 



which is exactly the Bohr-Sommerfeld quantization rule 
obtained recently in [ll|. The residues at each pole give 
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FIG. 2: Exact (top) and semiclassical (bottom) Husimi distri- 
butions for the 200th state of the LMG model, with j = 200, 
u> = 1 and Ua = 1000. 



FIG. 3: Exact (top) and semiclassical (bottom) Husimi dis- 
tributions for the 270th state of the LMG model. The param- 
eters are the same as in the previous figure. 
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the Husimi functions 



n n {z) 



i 



j 

x exp 



\z\ [T{E n ) + (dX S K/d£)\ R 
1 (£„ - f (z) + A{z) f ' 



(29) 



Here we have multiplied by the coherent states normal- 
ization factor (1 + \z\ 2 )~ 2: > . These functions are our main 
result. They are in general not normalized, and have 
a strong resemblance with the canonical semiclassical 
Husimi functions presented in 0. 

As a first example, we calculate the Husimi function 
for the simple case H — hujJ z . The Hamilton equations 
can easily be solved and give z(t) — e~ lu)t Zi and z{t) = 
£ iuj(t-T)- ^ which implies \z\ = u\z\. The SK correction 
is also very simple, with A(z) — hu>/2 and dlsx/dE = 0. 
The final result is 



n m (r) 



1 



r (j(l-r) + (l + r)(m + l/2)) 2 



4j> 



where m goes from —j to j and r = \z\ 2 . When properly 
normalized, this approximates the exact distribution, 



r j+rn 



(j + m)l(J — m)\ (1 + r) 2 i 



(31) 



quite well for large values of j, as we can see in Fig 1. 

Now let us turn our attention to a less trivial sys- 
tem. In the authors have shown that the semiclas- 
sical quantization condition (j28(l works very well for the 



Lipkin-Meshkov-Glick (LMG) model 



H = hcuJ z + ah 2 [J 2 - J 2 } 



(32) 



already at moderate values of j. We now consider the 
accuracy of the semiclassical approximation l|29|l for its 
stationary states. For small values of a the results are 
very similar to the previous case, so we consider only 
j = 200, to = 1 and ha = 1000. In order to display 
the results, we use the canonical coordinates (q,p) given 
in l[7)l. in terms of which the phase space is compact, 
q 2 + p 2 < Ahj. We show the exact and the semiclassical 
(normalized) Husimi functions for two different states in 
Fig 2 and Fig 3. The agreement is excellent. 

Summarizing, we have obtained a semiclassical approx- 
imation for the phase space representation of stationary 
states of spin systems. This was done by investigating 
the semiclassical Green function in the vicinity of real 
periodic trajectories. The accuracy of the result was ver- 
ified by comparing it with exact calculations for the sim- 
ple case H = TiloJ z and for the Lipkin-Meshkov-Glick 
model. Husimi functions are known to be good tools to 
study quantum chaos |13| , and an extension of this theory 
to more degrees of freedom would be interesting in order 
to approach chaotic systems (a trace formula for chaotic 
spin systems was recently obtained |l4j but, as already 
noted, taking the trace obliterates the information about 
the residues). Work in this direction is in progress. 
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